Dipole emission and coherent transport in random media III. 
Emission from a real cavity in a continuous medium 
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This is the third of a series of papers devoted to develop a microscopical approach to the dipole 
emission process and its relation to coherent transport in random media. In this paper, we compute 
the power emitted by an induced dipole and the spontaneous decay rate of a Lorentzian-type dipole. 
In both cases, the emitter is placed at the center of a real cavity drilled in a continuous medium. 
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It is known that the spontaneous emission rate, F, in 
a dielectric medium depends on the interaction of the 
emitter with the environment [H. This is so because the 
surrounding medium determines the number of channels 
through which the excited particle can emit. That is, the 
local density of states (LDOS). For practical purposes the 
understanding of the spontaneous emission rate of a fluo- 
rescence particle is of great interest in biological imaging 
Closely related is the process of emission/reception 
by nano-antennas [§]. Also, inhibition of spontaneous 
emission is expected to occur in photonic band gap ma- 
terials [1] as a signature of localization. 

In a previous paper general analytical formulae 
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FIG. 1: (a) Cermet topology scenario. The emitter with po- 
larizability qo is surrounded by a disconnected net of scatter- 
ers of polarizability Qs and correlation length ^. (fe) Simply- 
connected-non-contractible topology in which the emitter is 
placed at the center of a cavity of radius R surrounded by a 
continuous medium of dielectric contrast e. 

were derived for LDOS and F as functions of both the 
electrical susceptibility of the medium and the geome- 
try of the embedding of the emitter in it. There, three 
different cases are addressed attending to the nature of 
the emitter. In the first case, the emitter is a polarizable 
dipole of transition amplitude /i at resonance. In the 
second case, the emitter is a non-resonant dipole of po- 
larizability a with induced dipole moment e^aEQ, where 
Eq is a fixed exciting field. In the third case, the emit- 
ter is a fluorescent particle seated on top of a polarizable 
host particle. Attending to the embedding of the emitter 
and the topology of the system emitter-host-medium, we 
differentiated two cases |5|. In the first one, the emit- 



ter is equivalent, replaces or is placed on top of a host 
scatterer. In such a case, by consistency, the topology of 
the medium is that of a disconnected manifold -known 
also as cermet topology- in the sense that the medium 
is made of spherical inclusions and the emitter takes the 
place of one of those inclusions. Attending to the em- 
bedding, the emitter is placed within a virtual cavity but 
for the case it replaces a host scatterer. The spatial cor- 
relation of the emitter with the surrounding medium is 
given by the same correlation function which enters the 
electrical susceptibility tensor, x- Spatial dispersion in x 
is inherent to that topology. This scenario was studied 
in Q with an emitter seated on top of a host scatterer. 
A relation between the complex refraction index of the 
host medium and the spontaneous decay rate was found. 

In this paper we will focus on the remaining topologi- 
cal case. The emitter is seated at the center of a spherical 
cavity drilled in a homogeneous medium. This medium is 
continuous as seen by the emitter provided that the cav- 
ity radius, R, is much larger than any correlation length 
between the constituents of the host medium, f . Thus, 
the cavity is real. On the other hand, if the wavelength 
of the emitted field, A, is that long that it cannot re- 
solve neither the microscopical structure of the medium 
nor the cavity, the medium is also continuous as seen by 
the propagating field. Consequently, the host medium is 
characterized by a scalar susceptibility x which does not 
present spatial dispersion. Topologically, this setup cor- 
responds to that of a simply-connected-non-contractible 
manifold with a unique hollow sphere -see Fig[T] Regard- 
ing the emitter nature, two scenarios will be addressed 
for their potential interest in biological imaging and sig- 
nal processing. In the first one, the emitter will be a 
polarizable particle excited by an external fixed field of 
frequency much lower than any of the resonance frequen- 
cies of the emitter. The physical quantity to compute in 
this case will be the power emission, . In the second 
case, the emitter will be a polarizable particle which de- 
cays at some resonance frequency. Dipole induction will 
be encoded in the emitter itself as it will be assumed to be 
of Lorentz-type. Thus, the net effect of the host medium 
will reflect on a variation in the spontaneous decay rate, 
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r", together with a shift in the resonance frequency with 
respect to those values in free space. 

In the first scenario, the formula for the power emitted 
by an induced dipole reads 
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where ujq is the frequency of the exciting field Eq , with 
ko ~ luq/c much lower than any internal resonance fre- 
quency and ao = 4770^1^^ is the electrostatic polariz- 
ability of the emitter in vacuum, being ee(wo) its dielec- 
tric contrast and a its radius. For our topological setup, 
the 7- factors 7j^ || are given by Q 
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with fc evaluated at ko and where the cavity factors G_l || 
are 



(a) 



-= G^" 



H=G, 



with cos = k-k'. In these equations, hc{k) is the Fourier 
transform of the correlation function which describes the 
isolation of the emitter within a real cavity. For a spher- 
ical cavity of radius i?, hc{r) = —Q{r — R). In the 
above equations G± and G|| are respectively the trans- 
verse and longitudinal components of the renormalized 
propagator of the coherent -macroscopic- electric field 
in the host medium. They read, G\(k) = ~^ ^^l'} , ^ , 

^ ' k-^e±{k}—k^ 

r^jTT, k being a unitary vector along the prop- 



G||(fc) = - 

agation direction and A(fc) = I — k ® k being the pro- 
jective tensor orthogonal to the propagation direction. 
In these expressions, ej^ ||(fc) — 1 -f X-L.||(^)- Because no 
spatial dispersion is assumed in the continuous medium, 
€±_ = e|| = e arc functions of k only. G'^j^i^\^ are the compo- 
nents of the propagator of the electric field in free space 
with bare wave number k. They are given by the above 
formulae, Gj^, G||, with ej^ = e|| = 1. 

In the second scenario, the emitter nature differs from 
that in [Q]. In the present case, the dipole moment of 
the emitter is not a linear combination of a fixed and 
an induced dipole. Rather than that, the emitter is po- 
larizable in its own and its polarizability function a{k) 
gets renormalized as a result of iterative self-polarization 
processes. 



FIG. 2: (a) Feynman's rules . (fe) Diagrammatic representa- 
tion of the transverse and longitudinal polarization factors, 
27_L and 7||, as given by Eqs. (|4l5p . 
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In the above equation, the 7-factors are those in 
Eqs.dHD. In addition, the term \aQk'^^{2^f^} has 
been introduced to account for the internal resonance 
of the emitter in vacuum. As shown in (3], it plays 
the role of a regulator of the intrinsic ultraviolet di- 
vergence in g'"^ . That is, 3fi{27^^} = -^^^ where ag 
is the corresponding real electrostatic polarizability and 
fcn is the resonance wave number in vacuum. Following 
[8|, by parametrizing Eq.® in a Lorentzian (L) form, 
UL^k) = aQ/c^g^[A:^g^-fc2-ir"A:3/(cfc^gJ]~\ we can iden- 
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tify the decay rate of the emitter in the host medium as 
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where fc^es is a real non-negative root of the equation 
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= ao(fco/fcres)^ is the renormahzed electrostatic po- 
larizability and Tq = caok^/GTr is the in- vacuum emission 
rate. As argued in fs*], consistency with Fermi's golden 

rule requires ao = -^j^^i fJ- being the transition ampli- 
tude between two atomic levels in vacuum. 

In view of Eqs. (|ll2l3l9p . we can address in parallel the 
computation of W" and for each scenario described 
above. The computation reduces to the evaluation of 
27±(fc) + 7||(fc) for k — ko and k = k^es respectively. 

Following 0, we will decompose F" and VF" in 
transverse and longitudinal components, F" = 2F" 
and = 2PF" -t- respectively. Neglecting longi- 
tudinal resonances as in 6], only the transverse compo- 
nents contain propagating modes, Fp and Wp respec- 
tively, which are proportional to the coherent intensity 
and so to 3{27]'(fc)}. It is convenient to identify their 
specific contribution to F" and W"' for observational pur- 
poses in imaging. Propagating modes are those associ- 
ated to complex poles of the propagators in the inte- 
grands of eq.Q. We can read their contribution from 
the diagrams in Figl2l 
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where the first term on the right hand side {r.h.s) of 
Eg. pTj) contains the contribution of bulk propagation 
-transverse components of diagrams (&i), (62), (^4) in 
FigH 

All the above integrals in Eqs.(|i H7TTT|) can be per- 
formed analytically at any order in kR. However, for 
the sake of consistency with the approximations consid- 
ered so far, we take the limit kR <C 1 and keep leading 
order terms, 
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It is our choice to normalize W" to the propagating 
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power in vacuum, VFq = 

VF" reads 



-nr^M — W3 la l-gyP- Thus, 
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where 7^p(ao) = |1 + ifc^ao[27j_(fco) 7||(^o)]| ^ is 
the renormalization factor due to self-polarization cy- 
cles. Written this way, we recognize in the first term 
of the r.h.s of Eq.p^ the usual bulk term corrected 
by Lorentz-Lorenz (LL) local field factors [l3|, F'^-'^ = 

Fo5ft{|^^^^ v^}. The second term there includes cor- 
rections of order > x^. It is also remarkable that if the 
empty-cavity Onsager-Bottcher (OB) [3, local field 

factors are used instead, F*^^ — ^ o^{(y-:^jp[^ v^} does 

agree with the two terms of Eq. (fT3)) up to order x^. The 
terms in E q.lTil) are associated to absorbtion in the host 
medium (6l. lll| while that in Ea.(fT3|) corresponds to ab- 
sorbtion in the emitter. The propagating emission reads 

Wo7^p(ao)3{27^(fco)} 

IFo7^p(ao)[5^{V^}5^{(^)'} 

1, 
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As shown in at leading order, the second term 
in Eq. p^ subtracts from the total F^-^ the non- 
propagating, non-absorptive longitudinal contribution. 

It is worth comparing this result with that obtained 
in There, provided that e follows Maxwell-Garnett 
(MG) formula ^9] and neglecting absorbtion, 

2. 
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Eq.lfTB]) and Eq. pT|) agree at leading order in x- However, 
in |6(| the cermet topology was considered and the emitter 
was placed at the site of one of the host scatterers. The 
coincidence of both approaches at leading order bases on 
the fact that the limit fco^ <C 1, with ^ being the radius of 
the exclusion volume between host scatterers, is implicit 
in Eq. p7|) . Therefore, no spatial dispersion enters x ~ 
hence, MG holds- and the radius of the virtual cavity 
i? = ^ also satisfies koR <^ 1 as in the present case. 

It is straightforward to write Wp in function of the 
effective refraction index fi and the extinction coefficient 
K = (2fco^ext)~^i with l^xt the extinction mean free path, 
using the identities 3?{x} = — 1 — and 9{x} — 2nK, 
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For the sake of completeness, we give also the values of 
the non-propagating emission in terms of n, k and ao, 
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where Eqs. (|19l20l2ip correspond to the power absorbed 
by the emitter, the transverse non-propagating power, 
Wf^, and the longitudinal emission, W^, respectively. 
Further on, by invoking causality, Kramers-Kronig sum 
rules relate n and lext and variations of can be ob- 
tained as a function of ujq • 

Let us introduce next the self-polarization effect. This 
is done by multiplying all the expressions above for W" 
by the renormalization factor TZp{aQ). Neglecting ab- 
sorbtion both in the emitter and in the host medium, we 
obtain at leading order in koR, 
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where Vr is the volume of the emitter cavity. In contrast, 
the OB renormalization factor due to self-polarization 
reads 
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Eq.imi) and Eq. 

which multiplies 



differ just in the local field factor 
f "Ph ^Se'" '^^ sach expression. With this 
result we finalize the computation of in the first sce- 
nario. 

In the following, we address the scenario where the 
emitter is of Lorentzian kind and compute Egs.fSl fTU)) . 
We first solve for the renormalized resonance wave vec- 
tor hfQs-i 
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At leading order in koR, the first two terms in Eq. ((24|) 
equal the expression found in Q for an interstitial emit- 
ter. The additional terms in Eq. (|25p are relevant for 
strongly absorptive media. That may be for instance 
the case of nano-antennas in metallic media Q. Addi- 
tional solutions to the one above for kres are only admis- 
sible beyond the small cavity limit and/or the continu- 
ous medium approximation. That is, for kR > 1 and/or 
k£^ > 1, Eq. fTn]) may contain several real solutions. That 
might give rise to a splitting of the in-vacuum resonant 



peak provided such an splitting is broader than the line 
width. The origin of the resonance shifts of Eas. (|24l25p is 
in the virtual photons which dress up the polarizability a 
through Eq. ([8|) . Those photons trace closed spatial loops 
as depicted by the diagrams in Fig. ((2). They amount to a 
classical continuum of sates not to be confused with real 
localized photons. The resultant frequency shift must be 
rather interpreted as the classical analog to the quantum 
Lamb shift jisl . [lit ■ On the other hand, the possibility of 
a resonance splitting would be subjected to consistency 
with our perturbative renormalization scheme and linear- 
response approximation. 

Finally, we give the full expression for F" at lowest 
order in koR according to Eq.Q, 
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Eqs. (|24j|26p can be compared with the results of 

In summary, following Q, we have found analytical 
expressions for the power emitted by an induced dipole 
W" -Eqs. (|13ll4ll5|22|) - and for the spontaneous emis- 
sion rate F" of a Lorentzian-type emitter -Eqs. ()9l26p - 
from a real cavity in a continuous medium. In the for- 
mer case we compute both the propagating and non- 
propagating power emission -Eq. p8| and Eas. (|19ll20)) 
respectively- in function of the complex index of refrac- 
tion. In the latter case, we also compute the shift in 
the resonance frequency kres -Eqs. p0|24|25|) . We found 
that, in the small cavity limit, kR <^ 1, the tuning of 
kres may be sensitively affected by a strongly absorptive 
metallic environment -Eq. (j25| ). Beyond that limit and/or 
the continuous medium approximation, we argue on the 
possibility of obtaining additional resonance frequencies. 
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